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The W— fusion scattering process W + W~ — > ZZ for off-shell W bosons is studied, focusing on the 
issue of its high-energy behaviour which is known to be anomalous. It is shown that the unitarity 
violating terms can be isolated and extracted in a well-defined and efficient way using the pinch- 
technique. This restores the good high energy behaviour of the cross section and, in particular, 
makes possible the identification of the Higgs resonance in the invariant mass distribution mzz 
of the Z pair. The discarded terms, which are proportional to the off-shellness of the W bosons, 
^ ■ cancel against similar terms originating from the remaining diagrams for the full physical process 
fif 2 — > f\f' 2 ZZ. This cancellation ensures the gauge invariance of our result, which therefore 
constitutes a meaningful separation between signal and background when they both contribute 
coherently. Equipped with this result, we are able to define a resonant approximation for the 
process pp — > ZZ + 2 jets + X, which circumvents the problem of good high energy behaviour 
without having to resort to the lengthy calculation of the complete set of diagrams. In this 
approximation only the W— and Z— fusion signal graphs are included, i.e. the ones which contain 
the Higgs resonance. We have verified that the approximate resonant cross section describes very 
well the full result not only close to the resonance but also beyond it. 



1 Introduction 



The discovery of the Higgs boson is the primary physics goal of the LHC pp collider. For low 
Higgs masses the dominant production mechanism is gluon fusion, gg — > H, but for a heavy Higgs 
there is also a significant contribution from the W — fusion process, WW — > H, where the Ws are 
emitted from incoming quarks.[] Since a heavy (Standard Model) Higgs is most easily detected 
in the ZZ decay channel, the process of interest is WW — > H — > ZZ. However the s— channel 
Higgs resonant diagram is only one of several diagrams contributing to this scattering process, the 
other diagrams containing, for example, trilinear and quadrilinear gauge boson vertices. In fact 
the W — fusion process 

W( Pl )W(p 2 ) -> Z(h)Z(k 2 ) (1.1) 

provides a classic illustration of the subtle gauge cancellations encountered in non-Abelian gauge 
theories. The role of the Higgs graph is crucial in obtaining a cross section which is well behaved 
at high energies, see for example the discussion in 0. In practice, however, one actually has to 
consider the case of off-shell W — fusion, since V^s emitted from the incoming quarks have q 2 < 0. 
But when the Ws are off-shell the delicate gauge cancellations responsible for the good high 
energy behaviour of the amplitude are spoiled. Unitarity is badly violated by terms which are 
proportional to the off-shellness of the Ws [|J , 

W 1 = p\-M 2 w 

W 2 = pl-Ml . (1.2) 

Unitarity is only restored for the physical cross section when the full set of diagrams for the 
qq — > qqZZ process is taken into account 0,0, i.e. not just the subset containing WW — > ZZ. 
This involves a very large number of additional Feynman graphs. 

Since in the region of Mzz ~ Mr the on-shell WW — > ZZ proces is clearly well approximated 
by the s-channel Higgs resonance graph (at least for M H not too large so that Y H <C M H ), it is 
interesting to ask whether the full qq — *> qqZZ process, including off-shell WW scattering, can be 
similarly approximated, while retaining good high-energy behaviour. In other words, one seeks a 
minimal set of diagrams that include the Higgs resonance while at the same time not spoiling the 
delicate gauge cancellations. 

To begin with one could consider only the Higgs resonant diagram contributing to qq — > qqZZ. 
However this leads to a cross section that behaves badly at high energies, growing like s. As a 
result the shape of the differential cross section around the resonance cannot be trusted. As 
a next step one could extend the set to contain all six of the W — fusion diagrams of process 
(|1 . 1| ) . As already mentioned, there are potential cancellations at high energy in this case but 
the cross section actually respects unitarity only when both Ws are on-shell. In the off-shell 
case when the six diagrams are embedded in qq — > qqZZ, see Fig. [|, this is not true anymore 
and in fact the situation gets even worse. The cross section now grows as s 2 with the leading 
terms always proportional to the off-shellness of the Ws, i.e. W\ and W 2 . For high enough 
energies the Higgs resonance gets completely swamped by these off-shell terms and the differential 
cross section exhibits no resonance structure. Only when the full set of graphs (Figs. 1 - 3) is 

^ee, for example, the review in JjJ. 
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considered does one obtain a cross section that behaves well at high energies. In this set large 
gauge cancellations take place and the resulting cross section exhibits the same behaviour as for 
on-shell Higgs production through the process qq — > qqH, namely it grows slowly as £n(s) for 
large s. 

The bad high-energy behaviour of individual diagrams can be attributed to the following two 
factors: i) the presence of the trilinear gauge boson vertices, and ii) the fact that the Z bosons can 
be longitudinally polarized. Both these factors bring extra momenta in the numerator of Feynman 
graphs. In fact the polarization vector of an energetic longitudinal Z boson behaves as 



As a result of the Mz factor in the denominator, the cross section can grow as powers of s/M§. 
However the momenta introduced by these factors into the numerator of the Feynman graphs are 
at the same time the momenta which make possible a 'communication' between graphs with seem- 
ingly different propagator structures. It is this communication that eventually leads to the gauge 
cancellations between different diagrams that restore the good high-energy behaviour. Obviously 
the cancellation of gauge-dependent terms between different graphs will be even more pronounced 
in a general gauge, due to the presence of the extra momenta in the longitudinal part of a 
gauge boson's propagator in this case. 

Any attempt to approximate Higgs production using only resonant production via WW fusion 
(or indeed any other multi-gauge-boson production process) cannot therefore succeed unless the 
above gauge cancellations are correctly taken into account. The recently developed tree level pinch 
technique Q , , provides exactly the right calculational framework for addressing this problem. It 
allows a rigorous definition of a gauge invariant sub-amplitude which can be used to approximate a 
full scattering amplitude. For example in [|J, the authors considered the process e + e~ — > W + W~ 
and were able to establish good high-energy behaviour for each individual square or interference 
term of the three contributing Feynman diagrams. 

In this paper we will apply similar methods in order to isolate the W — fusion part of the qq — > 
qqZ process in the form of a squared subamplitude that is both gauge independent and respects 
unitarity. This subamplitude squared can then be used to provide a well-defined approximation 
to Higgs production based only on W — fusion that, as we shall see, works particularly well in the 
region of the resonance and above. It should therefore provide a useful analysis tool for simulating 
Higgs production at, for example, the LHC pp collider, especially in a Monte Carlo context where 
the compactness of our expressions for the subamplitude avoids a time consuming calculation of 
the full amplitude. 

The paper is in essentially two parts. In the following section we discuss in some detail how to 
implement the pinch technique in the context of the W — fusion process. Having arrived at the final 
expressions for the approximate scattering amplitude squared, we then perform numerical studies 
to compare our results for various distributions with those obtained from the full calculation. Our 
conclusions are presented in a final section. 
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2 A unitarity respecting amplitude for off-shell W fusion 

In this section we will consider in detail the process 



which obviously can proceed via W— fusion. We will show how it is possible to isolate the 
W— fusion part of this process in a well defined way, and arrive at a squared amplitude that is 
gauge independent and respects unitarity. 

The complete set of Feynman diagrams for this process exhibits a plethora of different prop- 
agators, and we can use these to classify the various types of diagrams. First, we separate all 
graphs into two categories according to whether the initial quarks annihilate (S graphs) or not 
(T graphs). The T graphs are those that contain t-channel W propagators. There are in total 34 
T graphs (in the Feynman gauge) and we may further separate them according to their fermion 
propagator structure or, equivalently, according to the particles that emit the final-state Z bosons. 
Thus if we require both Zs to be emitted from a gauge boson, we obtain the usual six W— fusion 
graphs. These will be denoted collectively by Tww 



and are depicted in Fig. 1. Their characteristic feature is that they do not contain any fermion 
propagators. In fact their structure is described by 



u(qi)u(q 2 ) -> d(q 3 )d(q 4 )Z(k 1 )Z(k 2 ) 



(2.1) 



Tww — Tw + T c w + T 4 g + T ( p + T c ^ + T# 



(2.2) 





(2.4) 



(2.5) 




(2.6) 



(2.7) 



(2.8) 




(2.11) 



(2.10) 



(2.9) 
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where an overall factor of (—ig 2 ) has been omitted from all graphs. 
The kinematics of the process is described in Fig. 1, i.e. 

Pi = qi ~ <?3 , P2 = <?2 - <?4 , Pi+ P2 = h + h = q ■ (2.12) 

The inverse bosonic propagators appearing in these graphs are 

Di = {k l -p l ) 2 -M 2 w = {k 2 -p 2 ) 2 -M 2 W) (2.13) 
D 2 = (h-p 2 ) 2 -M 2 w = {k 2 - Vl ) 2 -Ml , (2.14) 
D H = q 2 -M 2 H . (2.15) 

Note that an imaginary part must be included in the Higgs propagator Dh in order to regulate it 
when q 2 = Mfj. We will comment on its precise form later. The t-channel propagators of the Ws 
will always be spacelike (< 0). In fact for the two Ws emitted from the external fermions we have 

Pi = -^(1 - cosM , p\ = -^(1 - cosM , (2.16) 

where the energies and scattering angles refer to the centre-of-mass frame. Thus the Ws will always 
be off-shellby at least an amount equal to their mass, \Wi\, \ W 2 \ > M^, with the minimum being 
attained for forward scattering, #i 3 = = 6> 2 4 (=>- p\ — = pi). 

The forms of the trilinear Z^W~Wp and quadrilinear Z^Z„W~Wp gauge boson vertices that 
appear in these graphs are given, respectively, by 

rV*/?(?,Pl,P2) = (q-Pl)f39^a + (Pl-P2)i^9al3 + (P2-q)a9f3ix , (2.17) 

G a ^uipi,P2, h, k 2 ) = 2g a/3 g^ u - g a ^gp u - g au g^ , (2.18) 

where all momenta are considered incoming. 

The case where one Z is emitted from a fermion line while the other is emitted from a W is 
described by the eight graphs denoted by T W f 

Twf = T ld + T lu + T 2d + T 2u + T 1( j + T la + T 2d - + T 2n , (2. 19) 

and depicted in Fig. 2. These graphs contain one fermion propagator and one trilinear gauge 
boson vertex. Only one of the currents V a or Up appears. A typical Twf graph has the following 
structure: 

Ti d = t£- — r aw (pi, fci - pi, -k 1 )£ d ^=u d -f u P L - — —Yu u e^e 2 . (2.20) 

The couplings of the fermions to the Z are defined here as 

—-f(tfPL + r f P R ) (2.21) 

where 

if = J wj ~ slQf, r f = -s 2 w Q f . (2.22) 
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In the above expressions Iyyt is the third component of the weak isospin, Qf the charge of fermion 
/,andiV = (l± 75 )/2. 

The remainder of the T graphs, where both Zs emanate from fermions, are denoted Tff and 
are exhibited in Fig. 3: 

Tff = T Wl + T W2 + T Dl + Td 2 . (2.23) 



Each one of these graphs contains two fermion propagators. In Eq. ( |2.23| ) we have grouped them 



according to the W propagator they contain. In this class of graphs either one current appears 
and two fermionic propagators are present in the other fermion line (T^ and T\y 2 graphs, 6 graphs 
each), or no current appears and each fermion line contains one fermion propagator (Td 1 and Tjj 2 
graphs, 4 graphs each). Altogether there are 20 graphs in this class. 

Finally, the S graphs are divided into classes according to the neutral gauge boson (7, Z, or 
g) into which the initial quarks anihilate. There are 63 graphs in this category, denoted as Sz 
(20), if> 7 (20), S g (20), and Szz (3) according to the gauge boson propagators that they contain. 
In the Sv graphs the Zs are emitted from fermions and thus all these graphs contain two fermion 
propagators. The Szz graphs are the Higgs-strahlung graphs. Like the W fusion graphs they 
contain no fermion propagators. One representative of each of these classes is shown in Fig. 4. 

These seemingly different classes of graphs contain in fact a large number of identical terms. 
These common terms are responsible for the large unitarity cancellations that take place in the 
sum of all graphs. On the other hand, unitarity violation occurs within individual subsets of 
graphs, e.g. Tww> as a result of incomplete cancellations among such terms. The common terms 
among the different classes of graphs arise when momentum factors in the graphs' numerators 
trigger one of the following tree-level Ward identities: 

fc" TAt = # = (It-4i)+4 i = S- 1 {k-qd-S- 1 (q i ) 7 (2.24) 

q fl T lial3 (q,p 1 ,p 2 ) = \plg a /3 ~ P2aP2?} - [p\g a p - PlaPlp] = Uj(p 2 ) W ~ (pi)w, ( 2 -25) 

(pi,P2;h,k 2 ) = r a/3l/ (pi,p 2 + fa, fa) - r aj9l/ (pi + h,p 2 , fa), (2.26) 

& = \(q 2 -kl-kl)= l -D H + l -(M 2 H -2M 2 z ), (2.27) 

where Eq. ( |2.27|) is the Ward identity of the HZZ vertex with the two Z bosons on shell, while 
Uap(p)w — (p 2 — ^w)9a(3 — PaPp is the inverse propagator of the W in the unitary gauge. The 
inverse propagators generated in this way cancel one of the propagators of the graph, resulting 
in a structure that mimics the structure of a different class. For example, the two classes of 
graphs WW and Wf become identical in form when one of the Wi propagators is removed in the 
Tww graphs while at the same time the fermion propagator is removed in the T\yf graphs. The 
remaining terms in the Ward identities mostly cancel when the particles involved are on shell. 
As we have already noted, the momentum factors in the numerators that trigger the above Ward 
identities are furnished by the trilinear gauge boson vertices and the longitudinal polarization 
vectors of the external Z bosons. In fact current conservation, good high energy behaviour, the 
equivalence theorem and the cancellation of gauge parameters all have their origins in these and 
other similar Ward identities satisfied by the tree-level vertices of the Standard Model. Indeed a 
set of tree-level graphs is characterised as gauge invariant when none of its parts can resemble the 
structure of a different set of graphs by use of a Ward identity. 
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With this criterion in mind it should now be clear that the S graphs are not involved in any 
unitarity cancellations with the T graphs. First of all the couplings involved are different. For 
example the S g graphs can never mix with any other class of graphs since they are the only ones 
that involve the strong coupling constant. Thus they are gauge invariant and well behaved by 
themselves. Second, there are not enough factors in the numerators of these two classes that are 
capable of cancelling the necessary propagators to make a common term. One W propagator will 
always survive in the T graphs while no such propagator is present in the S graphs. Thus in the 
search for a well behaved W— fusion squared amplitude the S graphs can be ignored, and we need 
only concentrate on the T graphs. 

With these considerations it should also be evident that off-shell Z— fusion or Higgs-strahlung 
cannot give rise to unitarity violation ||. The Higgs-strahlung graphs Szz cannot communicate 
in any way with the rest of the graphs. Even if their Higgs propagator is cancelled they will always 
contain two Z boson propagators. 

Equivalently the good high energy behaviour of the T graphs alone can also be verified by 
considering a different W^-fusion process where the initial fermions cannot annihilate, like e + — > 
vv^ZZ . This process still retains the T graphs but does not involve any S graphs. 



The Ward identities of Eqs. (|2.24J - |2.27D guarantee both the gauge parameter independence of 



the amplitude as well as independence with respect to gauge transformations of the external gauge 
bosons. For physical amplitudes with external gauge bosons, gauge invariance is encoded in the 
following relationships: 

T ,g: WT^-g) = (2.28) 
W ± : PT At (H/ ± ) = ±M w T((f)) (2.29) 
Z : fe"T M (Z) = -iM z T{ X ) (2.30) 

In these equations 7^(7; g), T /J (W /± ), T M (Z) are physical amplitudes with at least one external 
gauge boson 7, g, W ± , Z carrying momentum k^. The amplitudes T(<f>) and T(x) 011 the right- 
hand side are the corresponding set of graphs where the gauge boson is replaced by its would-be 
Goldstone boson. The above Eqs. (|2~28| - ^30D can be considered as Ward identities for whole set of 
graphs and they are the direct consequence of the Ward identities of Eqs.( |2.23| - [2.27| )A set of graphs 



that satisfies these Ward identities is said to be gauge invariant. It would be both independent of 
any gauge parameters and of the choice of polarization vector for the external gauge boson. For 
the photon, Eq. (]2.28|) is usually referred to as current conservation in QED, while from Eqs. (|2.29| ) 



and ( |2.30| ) the Equivalence Theorem follows directly [|J. 

In our case, since the fermions are massless and there is no \WW vertex, the T\yf and Tff 
graphs will give no contribution to the right-hand side of Eq. ( |2.30| ). Thus 

KT^Z) = -iM z T ww ( X ) . (2.31) 

Since unitarity in off-shell W— fusion is violated by terms which are proportional to the off- 
shellness of the W bosons [0, the common terms between the T ww and the T W f + Tff graphs will 
always be proportional to W\ and Wi- In the calculation of the squared amplitude such common 
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terms will emerge at three different levels. The first level is that of the amplitude, stripped of the 
polarization vectors of the Zs 



T>» = (T ww + T Wf + T ff Y u . (2.32) 

At this level, the presence of the trilinear vertices already generates off-shell unitarity violating 
terms in the Tww and Twf that cancel between them. At this stage all gauge parameters cancel 
also. 

The second cancellation occurs when in the full amplitude, 

T = T^e^ , (2.33) 

the polarization vectors of the Z bosons become longitudinal and thus their leading part becomes 
proportional to the momenta of the Zs, fcf and k 2 . Now the action of fcf and/or k 2 on the 
amplitude will produce extra off-shell terms. Such terms will again cancel among the Tww and 
T W f + Tff graphs. 

Finally, in the squared amplitude additional off-shell terms will appear between the squared 
W— fusion amplitude iTy^! 2 and the interference term T WW {T W ^ + ^7/)- Most of these terms 
again cancel. Those that do not, emerge from the interferences of the Higgs graph T#. Because 
of their structure (they contain no fermion propagators), these terms will be allocated to the 
W— fusion amplitude |7\yjy| 2 . 

Accounting for these cancellations at every level is simply an exercise in identifing structures. 
The common terms can actually be identified graphically, as in Fig. 5 or Fig. 6. In order to keep 
track of them the properties and the Ward identities of the gauge boson self-interaction vertices 
are exploited. Thus the trilinear gauge vertex of Eq. (|2.17| ) is decomposed in the asymmetric form, 
first used by t' Hooft [|], 

r Ma/ 3(g,Pi,p 2 ) = r£ a/3 (g;pi,p 2 ) + T^{q; Pl ,p 2 ) , (2.34) 

where 

r JL/3(<?; P11P2) = (pi - P2)^g a p - 2qagp» + ^g m (2.35) 

T l a p((liPl,P2) = -Plagp»+P20g m ■ (2.36) 

The pinch part is proportional to the momenta of two of the legs of the vertex carrying their 
corresponding Lorentz index. This part therefore vanishes identically when the two legs are on 
shell. The F F part of the vertex satisfies the Feynman gauge part of the Ward identity of Eq. (|2.25|) , 

q^(q; Pl ,p 2 ) = (pj-p^g^ = (W 2 - W 1 )g a p , (2.37) 

while the pinch part T p will give rise to the longitudinal terms of Eq. ( |2.25| ). The four-gauge-boson 
vertex of Eq. (|2.18f) can also be split into two parts: 

G a p-,nv{pi,P2, -ki,-k 2 ) = G^ u (p 1 ,p 2 ;-ki,-k 2 ) + G^ v (p l ,p 2 ;-k l ,-k 2 ) , (2.38) 
G^ v {Pi,P2;-h,-k 2 ) = (2.39) 
G^M,P2-,-ki,-k 2 ) = -sTg^-sT^. (2.40) 
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With the above decomposition, the Ward identity of Eq. ( |2.26|) is also satisfied individually by its 
Feynman (G F ) and pinch (G p ) parts. 



We next proceed to explicitly exhibit and implement the aforementioned cancellations between 
the W— fusion and the rest of the T graphs. Most of the remaining part of this section will be 
quite technical, involving extended use of the relevant Ward identities. Readers who are only 
interested in the results of the calculation may wish to go directly to Eq. ( |2.92| ) and the discussion 
that follows it. 

The 1st cancellation 

In order to identify the common off-shell terms at this level one must first decompose the 
trilinear and quadrilinear gauge vertices into their Feynman and pinch parts. The two sets of 
graphs Tww and Ty/f which contain such vertices will then separate into Feynman and pinch 
parts respectively: 

rpflU rpFflU , rpPplU 

1 ww — 1 WW ' 1 WW j 

TVs = T w f + T^/. (2.41) 

The pinch parts are proportional to the off-shellness of the W bosons. We will explicitly show 
that these parts exactly cancel. 

Since the decomposition of a trilinear vertex in Eq. Q2.34Q is asymmetric, at this point one 
has to make a choice for the first momentum argument of the T F vertices. Different choices 
should make no difference to the final result when the procedure is carried through consistently. 
Nevertheless at intermediate points a particular choice may facilitate the calculations. It turns out 
that choosing as the special momentum for the T F vertices the momentum of the corresponding 
off-shell W, pi a or p 2 p, minimizes the number of interference terms that must be considered, thus 
expediting the calculations.^ Thus the Feynman and pinch parts of the vertices are: 



T app(Pu fa-pi, -fa) = 2(ki a g pli +pipg ap -Pxp9iux) ; r^Oi; fa- 

t„ .1. „ u \ _ oa „ , „ „ „ „ \ . rP 

Pp, 
ap 

^F /„ . ;„ „ u \ _ nlu „ i „ „ „ „ \. rP 



r p P v (P2 ; h - P2 , - h ) = 2 ( ktpgpy + P2v9pp ~ P2 P g v p ) ] V J pjy ^ 2 , ^2 ■ 
FapviPu h-pi,-k 2 ) = 2(k 2a gpv+pivg a p-pipgu a ) ; rj»^(pi; h- 



Pi, 


-h) 


= (Pi 


kl) p gpa 


(2.42) 


P2, 


-k 2 ) 


= (P2 


— k 2 ) p g U {3 


(2.43) 


Pi, 


-k 2 ) 


= (Pi 


— k2) p gva 


(2.44) 


P2, 


-hi) 


= (P2 


-ki) p gpi3 


(2.45) 



T /3pp(P2; h-p 2 , -fa) = 2(ki f3 g pil +p 2ll g l3p -p 2p g^) ; Y pp ^ 

In the above expressions we have dropped all terms that will not contribute to the matrix ele- 
ment. These are the terms proportional to pi a and p 2 p, which for massless fermions vanish when 
contracted with the external fermionic currents, i.e. pi a V a = and p 2 $U^ = 0. Furthermore at 
this level we can also drop terms proportional to fa p and k 2v that vanish when contracted with 
the polarization vectors of the Zs: ki p e^{ki) = = k 2v e u {k 2 ). 

It is straightforward to calculate the pinch part of the TV/ graphs. Since they contain only 
one trilinear vertex they split immediately into two parts. One readily observes that the pinch 

2 Actually this choice is identical to the one made in the one-loop pinch technique, since these would be the 
momenta outside the loop in WW — » ZZ. 
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part, T p , of the vertex simply cancels the fermion propagator in each of these graphs by virtue of 
the Ward identity of Eq. fl2.24|) . We illustrate this explicitly for one of the graphs. The expression 
for Ti d of Fig. 2 has already been written out explicitly in Eq. fl2.20|) . Using Eq. ( [2.42|) for Tp , 
the pinch part of this graph is given by 



Then writing 



T[ d = wTT^^a T^T^ ~ h) p L««(fb)$<* ■ (2-46) 

W\ Di V2 ii -Ipi+fi 



(2.47) 



the first term will cancel the fermion propagator while the second one vanishes since ^uj^li) = 0. 
So the pinch part of this graph is 

V a U p Wo 

T " = wiw^ir^^ ■ (2 ' 48) 

The corresponding expression T^, of graph T\ u of Fig. 2 is simply obtained by (i) changing £ d , 
the fermion coupling of the Z, to l u since the Z boson is now emitted from the up fermion, and 
(ii) changing the overall sign, since now the fermion propagator involves ^3 — ($ii — ^1). Thus, 
using the following relation for the fermion couplings, 



ct , (2.49) 



the two graphs combined give 



V a Up 2 .W 2 



Tf = T[ u + T[ d = ^^(-o^ga.gpu ■ (2.50) 

Evidently this expression resembles the structure of the W— fusion graphs, T^w- It contains no 
fermion propagators and has the correct coupling. We have divided and multiplied by the W 
boson inverse propagator W2 in order to bring the expression into the form of Eq. ( |2.3|) . This step 
is not strictly necessary, but our convention will be always to extract a factor V a Up/ '(W1W2) from 
all such terms. 

The pinch parts of the rest of the Ty/f graphs are extracted in a similar way. The graphs 
always combine in pairs to produce the correct coupling (? w . Altogether we obtain for the pinch 
terms of the TV/ graphs the following expression: 

= wM { ~ Cl){Wl + W ^ 9 -^ + ^K £ 2 • (2-51) 



Next we turn to the W— fusion graphs. In order to identify similar pinch terms to those of the 
T W j graphs we will use the following identity for the product of two trilinear vertices 0: 

T api jTpp U = T app Tp pu + T apfJ Tpp V + T app T p pu — T app T p pv . (2.52) 
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This decomposition will enable us to make use of the Ward identity of the full vertex given in 
Eq. (pi). 

Let us first consider the graph T w given in Eq. (|2.6|) . Using Eq. (|2.52|) , the product of its 
trilinear vertices can be expressed as 

r apAt (Pi, h-p 1} -k^Tppvfo, k 2 - p 2 , ~k 2 ) = rJ, M (pi; h-pi, -k 1 )Y^ pv {p 2 ; k 2 - p 2 , -k 2 ) 

+ gpu(P2-k 2 ) p r a ^( Pl ,k 1 - Pl ,-k 1 ) 
+ 9a,x{pi - fci) P r^(p 2 , k 2 - p 2 , -k 2 ) 

- 9aix9py{P2 - k 2 ) ■ ( Pl - ki) . (2.53) 

Since p± — ki = k 2 —p 2 , the Ward identity of Eq. (|2.25|) is immediately triggered by the second and 
third term of the above equation. This will produce off-shell terms proportional top 2 = Wi + M w . 
The fourth term can be written in terms of the inverse propagator D\ when M w is added and 
subsequently subtracted. Finally, using the on-shell condition for the Zs, k\ = k\ = M|, the 
expression ( |2.6| ) for the diagram Tw is transformed into 

" + i M z( c l - 2s D + °i + (Wk + W 2 )}g a »gJ . (2.54) 



I 



w D 1 

The momenta arguments of the vertices will no longer be exhibited, since the way their momenta 
are assigned should be obvious from the Lorentz indices of the vertices. The first and third 
indices determine the first and third momentum arguments of the vertex respectively, according 
to (/!, u, a, (3) — > (—hi, —k 2 ,pi,p 2 ), while the second argument is fixed by momentum conservation. 
In an identical way the crossed graph T c w is rewritten as 



2 



rpOlfijiV Cy) 



rL,rL, + [m 2 z {c 2 w - 2 S 2 W ) + d 2 + {w x + w 2 )} gau9 ^ (2.55) 



Having obtained these new forms for the diagrams Tw an d T c w we make the following ob- 
servations. The terms proportional to M|(c^ — 2s^), directly combine with the corresponding 
would-be Goldstone graphs of Eq. (|2.9|) and Eq. (|2.10| ) to produce an overall coupling equal to 
s w ~ 2s w c w + c w = ( c w ~ s w) 2 - The terms proportional to D\ or D 2 will immediately cancel 
the relevant propagator of the graph, i.e. either 1/Di or 1/D 2 , and will thus combine with the 
quadrilinear vertex graph T^q. In doing so they cancel the G p pinch part of the quadrilinear 
vertex. The last terms proportional to the off-shellness of the W bosons are the pinch terms. 
These operations are represented pictorially for graph Tw in Fig. |](a), where the ellipsis represent 
the M 2 z \cl - 2s 2 J term of Eq. ( ggg) 

Thus the Feynman and pinch parts of the Tww graphs are given by the following expressions: 

V TT r 2 r 2 
rpFiiv _ v a Up f C w F p C w p p 2 

1 WW — yy^ yy^ \ £) a PP PpV £) 2 apU ^w9aP9ilV 

+ M 2 z {c 2 w - s 2 w ) 2 p^ + ^SH] + M|^^| , (2.56) 
D\ D 2 D H > 



and 



-,Pfiu _ Va Up, 2 ,, \r9an9Pu 9av9p^ 



T7w = wM { J{Wl + W ^^BT + ■ (2 ' 57) 



10 



We note that the off-shell pinch terms have only been generated from the graphs that contain 
trilinear gauge vertices, namely T w and T cW . Only these two graphs can communicate with the 
rest of the graphs, in this case Twf, due to the rich momentum structure of their numerator. 
Finally we observe that the pinch part of the W— fusion graphs, T^/w m the above equation, 



exactly cancels the pinch term of the TV f graphs T w j of Eq. ( j2.51| ) . Indeed 



T W w + T W f = . (2.58) 

We have also verified the above results by performing the calculations in a general gauge. 
For this 2 —>■ 4 process the cancellation of the gauge parameters due to current conservation is 
not automatic as in the case of a 2 — > 2 process. Common pinch terms among different classes of 
graphs are now even more prolific, due to the presence of the extra momenta in the longitudinal 
part of the gauge bosons' propagators. In addition, most of them will also be gauge parameter 
dependent. When all such terms are identified and cancelled the surviving expressions coincide 
with those obtained in the Feynman gauge. Thus the formula of Eq. ( |2.56| ) for T ww and the 
corresponding ones for T W j and Tff are truly gauge parameter independent expressions. 

It would be interesting to find or invent a gauge where these expressions could be obtained 
automatically by the Feynman rules of the particular gauge. In such a gauge this first level of 
cancellations would be avoided. Furthermore such a gauge might prove useful in other multi- 
gauge-boson processes. However we are not aware of any such special gauge. Although the T F 
vertices of Eqs. ( p.42| - [2T45| ) look similar to the trilinear vertices of the Gervais-Neveu gauge f9j, 



they are actually not the same, and thus the Gervais-Neveu gauge expression of the W — fusion 
graphs does not coincide with T ww . 

Although the expression for the W fusion graphs Eq. ( [2.56D , obtained after this first step, is 
gauge parameter independent, it remains however gauge non-invariant in the sense that it still 
does not satisfy the Ward identity of Eq. ( |2.30| ). Because of this fact, the W fusion part can still 
not be separated from the rest of the graphs. To do this, the cancellations inherent in Eq. ( |2.30| ) 



must be allowed to take place first. This is done in the next step. We also note that the following 
steps would not be necessary had we considered photons instead of Zs in the final state, i.e. 
WW -> 77. 

The 2nd cancellation 

After this first cancellation has taken place the resulting gauge parameter independent ampli- 
tude T for our process can be written as 

^ = C^ww + T ww + T W j + T Wf + TffY u eie u 2 

= (T ww + T^ f + T ff re^ 2 . (2.59) 

The further cancellations due to the longitudinal Z bosons become more transparent in the 
squared unpolarized amplitude where the sum over the polarizations of the Zs will give a factor 

£4(*i,A a )4" (ki, A x ) = -g^ + ^-f , (2.60) 

Ai 1V1 Z 
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and similarly for e z (k 2 , X 2 )- The extra momentum factors, k%, k% etc. introduced in this way will 
result in bad high energy behaviour, since their growth with energy cannot be compensated by the 
constant factors of M z in the denominator. Thus cancellations at this stage will be instrumental 
in restoring unitarity. 

The unpolarized squared amplitude is given by 

~W= -T F ^(-g, p , + A^Av/Mf ) (-g w , + k 2v k 2v , / 'M 2 z )t f *, . (2.61) 
Expanding the product of polarization tensors gives four terms: 

1 1 ^ M z M z 



M z ' Mz'^ Ml M 2 Z 

We next determine the effect of the factors of longitudinal momenta within each of the following 
terms: 

k lp T F ^ } k 2v T F »\ and k 2u k lp T F ^ } (2.63) 

before actually squaring them. We will then explicitly show the generation and cancellation of 
the common off-shell terms within each of the above terms. 

The action of k^ on the T ww amplitudes will generate off-shell terms when contracted with 
the T F vertices. Since k^ is not the first, special, argument of T F the Feynman Ward identity of 
Eq. ( ^.37| ) is modified to: 

K r£ PM = (~Di + Mf + W x ) g ap + 2k la (h - Vi) P (2.64) 
for the graph 2V, and to 

k^ F pp = (-D 2 + M 2 Z + W 2 ) g Pp + 2k ip {k x - p 2 ) p (2.65) 

for the crossed graph T c w- The remaining terms (k\ — p\) P and {k\ —p 2 ) p of the above equations 
will create additional off-shell terms when they act in turn on the remaining T F vertex of each 
graph. Now the modified Ward identities read: 

(A* - Pl yr F pu = — [Di + W 2 + Ml(4 - sl)]g Pv - 2k w k 2v , 

{h - p 2 ) p r F p „ = _[£> 2 + Wi + M\{cl - sl)]g au - 2k 2a k 2v , (2.66) 

where the following elementary identities have been used for the dot products: 

2h ■ Pl = -D x + M 2 z + W l , 2kt ■ p 2 = -D 2 + M 2 z + W 2l (2.67) 
2A; 2 • Pl = -D 2 + M 2 Z + W 1 , 2k 2 ■ p 2 = -D x + M 2 Z + W 2 . (2.68) 

In the end the contraction of /cf with the T^r W graphs assumes the following form: 

kOTX = M l A i - 24^i*a + M%^-k\ + c 2 w O\ . (2.69) 
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In the term A\ the index v is never carried by either hi or k\. The dependence on k% or k\ has 
been explicitly exhibited, the latter emerging from the Higgs graph. The explicit expressions read: 



V a U p 

w[W 2 

v a u p 
W[W 2 

v a u 13 



A\ 
Ft 



t «/3 



D 2 



F 

a(3v 



(2c 



r 



D\ ^ ~D 2 



W x W 2 



a 3 



with 



,af3 



a(3 , Zfi a K 2 



CJ + 



+ 



Do 



The off-shell terms 0\ in (|2.69|) are 



V a U p 

= 1 

1 W x W 2 



Wx P Wo „ Wo Wx ft ' 



£>9 



Dx 



Do 



(2.70) 
(2.71) 
(2.72) 

(2.73) 
(2.74) 



The action of k± on the TV/ and Tff graphs can be easily determined. Since mostly the Z bosons 
are emitted from fermion lines, the k± will directly pinch the fermion propagators adjacent to the 
Z boson leg by virtue of the Ward identity of Eq. ( |2.24j ) with k — ► kx- 

When ki is emitted from a W then the modified Ward identities of Eqs. Q2.64|) and (|2.65| ) 
must be used. The terms (kx — Px)p and (kx — p 2 ) P emerging from these identities will now pinch 
the fermion propagator attached to the internal W boson the same way as in Eq. (|2.46|) . From 
the remaining terms of equations (|2.64|) and (|2.65|) the ones proportional to M§ cannot obviously 
pinch while the rest pinch one of the W propagators of the graphs. These latter pinch terms attain 
therefore a structure similar to that of the Tff graphs and indeed will cancel exactly the entire 
contribution of the Tff graphs. 

After carrying out the above steps we finally obtain 



K-l J- Wf 



-clO\ + M\R\ - klT 



f.f 



(2.75) 



One immediately observes that the off-shell pinch terms that emerge from the Tyyf graphs are 
exactly opposite to the ones generated by the W — fusion graphs of Eq. ( |2.69| ) and exactly cancel. 
This is a manifestation of the Ward identity of Eq. ( [2.30|) . The remainder R\ retains the structure 
of the Twf graphs and is given explictly by 



w 2 D 2 V2 [ 43 + h n ~ f 2 j 

Thus the action of k x on the complete amplitude will finally be of the form: 



r2 F ? 



kx^ = MIA\ - Ic^Fxkl + Ml-^ + M\R\ 



^2.76) 



(2.77) 
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As already mentioned, this is nothing else but the Ward identity of Eq. (|2.30|) slightly modified by 
the remainder of the Tyyf graphs. The fact that such a term has survived in the left-hand side is 
simply because after the first cancellation the amplitudes have been written in terms of modified 
vertices T F . Since we also do not contract with e u (k 2 ), the terms proportional to k 2 have also 
survived. 



The effect of the longitudinal factor k 2 in the third term of Eq. (|2.62|) is obtained in an identical 
manner. As expected, the off-shell terms common between the Tww and Twf graphs again cancel. 
In the end one obtains 



ko,T^ 



M 2 Z A% 



2^% + M 2 *k% + M Z R% 



D 



2 • 



H 



(2.78) 



where A 2 and R 2 can be obtained from Eqs. ( p. 70 ) and ( 2.761 ) respectively by the replacements 
ki <-> k 2 and fi <-» v. 

Finally we examine the fourth term of Eq. (p.62 ), k^ki^T^ . Here the momenta of both Zs 
act on the amplitude. (The order in which the contractions are performed does not matter.) 
Acting first with k± again produces the result of Eq. Q2.771) . When next k 2 is contracted with this 
equation it generates off-shell terms from both classes of graphs. On the Tww graphs it gives 



^2^1 1 WW ~ 1V1 Z 



, 2 4c v 



)F 1 + AF 2 + c 2 w 12 F 2 



with 



and 



A 



W 2 



Di 



+ 



M 2 H - 2M| 



2D 



H 



O 
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W l 

~D\ + D~ 2 



(2.79) 



(2.80) 



(2.81) 



Here we have also used the Ward identity of the Higgs vertex, Eq. ( p. 27f ), to extract a piece from 
the Higgs graph and combine it with the gauge boson graphs in F\. 

The result for the remainder of the Tjyj graphs R\ can be readily obtained by removing the 
fermionic propagators in Eq. (|2.76|) with ty 2 , by virtue of the Ward identity of Eq. (|2.24j) . This 
gives 

k\R\ = -ciO, 
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Finally collecting together Eqs. ( |2.79|) and ( p.82|) we obtain 



Notice again that all off-shell terms have cancelled. 



l-Acl)F 1 + F 2 A 



(2.82) 



^2.83) 



The 3rd cancellation 

After the first two cancellations have taken place, the unpolarized squared amplitude of 
Eq. ( p.62j ) attains the following form 



4 T 



rjlFflUrpF* 
1 )1V 
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M 2 



x 



M 2 

M\A\ - 2c 2 w F x k v 2 + F 2 -^k\ + M\R\ 

-2 

r2 av <-> 2 ui uv 1 17 ^"-2 1 n/f2 



M^* - 2<F 1 A£ + F 2 -^k\ + M|i^ 



D 



H 



M 2 . 



M 2 

( 1 --4c 2 w )F 1 + F 2 a\ 2 . 



(2.84) 



We can now see the advantages of our choice for the T F vertices. When each term of Eq. ( |2.84j ) 
is squared, additional off-shell terms will appear due to the interference of the T ww with the 
Tyyf + Tff graphs. This comes about from the action of the terms k x , k 2 , k x and k 2 which have 
been written explicitly in Eq. ( |2.84 ). No such terms will be generated in the first term however. 
Squaring the second term and exhibiting only terms that are of the W— fusion type we obtain 



1 



Mf 



M Z A X ■ A\ + 4ctMl|F 1 | 2 + Ml\F 2 /D H \ 



n v Ml u , , , ,2 



2ciF x k v 2 M Z A\ + M z F 2 -^k\ + M Z R\ + c.c. 



D 



H 



M 



+ F2- F rK { M l A \ + M 2 Z R V X + c.c. + 



(2.85) 



On the second line of this expression fl2.85|) , k 2 will generate the same off-shell terms, X2 , in a 
manner identical to that of the previous subsection. These terms will cancel in the same way they 
did in Eqs. ( |2.79|) and fl2.82j ). Thus the second line of (|2.85|) is equal to 



-2c 2 w F 1 M 2 z [(--2c 2 w )F* + A*F 2 * 



(2.86) 



Next we consider the third line of Eq. ( |2.85|) . This corresponds to the interference of the contracted 
Higgs graph k x T^ with whatever has remained. The index v is now carried by the four momentum 
k\ instead of k 2 . This will not affect the off-shell terms generated by the W fusion graphs which 
will be again X2 , 

M\k\A\ = M 2 z h, + M 2 z c 2 w 12 F 2 , 



(2.87) 



with 



V a U 13 

w[W 2 



D x D 2 



~ C w9af} 

+ (jgafiM^ + (1 - 4c 2 Jk la k^ 



+ 2c 2 w (ki a k 2p - k 2a k 



1/3) 



1 



1 

D~2 



{2.1 
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On the other hand, k\ acting on the remainder of the TV/ graphs R\ now generates off-shell terms 
with the opposite sign from the ones that would be generated from k 2 . This can easily be seen 
by looking at the expression for R\ in Eq. ( |2.76 ). For ^ to pinch the propagator 1/(^4 + ^2) we 
must use momentum conservation to write ^ = — k" 2 +$1 +$2, thus introducing a minus sign and 
a remainder proportional to j^i + ^2 = ^1 + ^2 which can give zero pinch. The final result will 
assume the form 

k\R\ = c 2 w M 2 z 12 F 2 . (2.89) 

Thus we observe that the off-shell terms in this case do not cancel between Eqs. ( j2.87| ) and ( |2.89| ) . 
The third line of ( |2.85|) is then equal to 



D 



H 



M 2 z hl + 2-c 2 w M 2 z 12 F* + .. 



(2.90) 



Since the off-shell terms that have survived in this interference of the Higgs graph do not contain 
any fermion propagator, they must be included in the expression for the squared W— fusion graphs 
and will therefore be a part of our final result. (Actually their numerical significance in the final 
result is negligible.) Similar results are also obtained for the third term of Eq. (|2.84|) with the 
replacement k\ <-> k 2 . 



Finally we collect together all terms that do not contain any fermion propagators, that is, 
all those that resemble the structure of a squared W— fusion graph, and omit all others. In this 
manner we are able to define a part of the squared matrix element specific to W— fusion: 



\T\ 



\T\ + 



(2.91) 



This squared amplitude will play the role of the signal for W— fusion, while all remaining squares 
and interferences will be identified as background. We stress again that this separation between 
signal and background is now meaningful since both parts are gauge invariant and well behaved 
at high energies. The signal squared amplitude for W— fusion is explicitly given by 



t f»v T *f _ M 2 Ai . A * 



M 2 7 A 2 ■ A% 



(1/4-84)1^1 



+ |F 2 | 2 [\A\ 2 - 2M Z /\D H \ 2 - AM^O/IDh] 



where h = hi + h 2 and 



+Ue{F 1 F%A*) - 2Mf 5?e {hF 2 /D* H ) 



O 



(2.92) 
(2.93) 



(q 2 - + W 2 )(l/D 1 + 1/D 2 ) . 

We observe that the O term changes sign above and below the resonance as an interference term 
usually does. Finally, the unpolarized W— fusion cross section can be calculated from 



9 A Vr\ 2 1 jk 



(2.94) 



Since in arriving at the result of Eq. ( p.92|) we have accounted for all possible exchange of 
terms and cancellations between W— fusion and T\yf +Tff graphs, this expression should now 
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respect unitarity. By direct numerical calculation (see below) we can see that this is indeed 



the case; Eq. (|2.94|) exhibits good high energy behaviour. In fact the cross section assumes a 
slow logarithmic growth with increasing scattering energy, exactly as in the production of an on- 
shell Higgs. Furthermore the Higgs resonance is clearly exhibited in the differential distribution 
da/dM zz where M zz = (£4 + k 2 ) 2 = q 2 . 

Since the kinematics allow q 2 to be equal to Mjj, in order to carry out the numerical evaluation 



of Eq. ( |2.94j ) the propagator of the Higgs boson must be regulated. We adopt here the naive 



approach of adding an imaginary part, equal to the total decay width of the Higgs, to all the 
Higgs propagator denominators in Eq. ( |2.92| ). Specifically, we make the substitution 

D H = q 2 - M 2 H + iM H T H . (2.95) 

The width T# is taken to be the sum of all partial, constant decay widths of the (Standard Model) 
Higgs boson, each one contributing as the relevant threshold is crossed. With all fermions apart 
from the top quark taken as massless, this is given byf\ 

T H = 6{q 2 - Am 2 ) rgP + % 2 - AM 2 W ) T^ W) + 6{q 2 - 4Mf) Y% Z) . (2.96) 

Including a constant width in this manner is to some extent an ad hoc assumption, but we have 
found that it works very well numerically. 

For a more consistent approach one has to rely on field theory. In field theory, regulators 
for resonant amplitudes are naturally provided through resummation of the self-energy diagrams 
which form a geometric series. The minimal approach, sufficient to regulate the resonance, is to 
include only the imaginary part of the one-loop self energy in the resummation. This corresponds 
to replacing the inverse Higgs propagator by 

D H = q 2 - M 2 H -> q 2 - M 2 H + i^mU HH (q 2 ) (2.97) 

thus giving rise to a running width 

^mU HH (q 2 ) = T H (q 2 ) . (2.98) 

Since the bosonic parts of the Higgs self energy are gauge dependent, the pinch technique must 
be used in both the construction and resummation of the self energies fL0| . This results in a 
gauge independent running width that coincides, at q 2 = Mjj, with the physical decay width of 
Eq. ( |2.96| ), term by term, viz. 

QmnP H (M%) = M H rP 1 -smIifP(M 2 H ) = M H Y^ w \ QmU^(M 2 H ) = M H T^ Z) . (2.99) 

The resummation of the Higgs self energy alone distorts the Higgs related Ward identities such 
as Eq. ( 2.27Q . This is because now the inverse Higgs propagator contains terms of 0(a). To 



compensate for this and still maintain the Ward identities, we must also include in our amplitude 



3 Recall that we are interested only in heavy Higgs bosons, decaying to Z°Z° final states, whose total width is 
dominated by the ti, W + W~ and Z°Z° contributions. 
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the one-loop imaginary parts of the relevant Higgs vertices. The full form of the Ward identity of 
Eq. flT27D is @ 



*$W% ZZ {q, ~h, ~k 2 ) + M 2 z T H ™(q, -fa, -k 2 ) = ^ [D H (q 2 ) - D x {k\) - D x (k 2 2 )] , (2-100) 
where x is the Goldstone boson related to the Z gauge boson. A similar Ward identity holds for 



the T^J VW vertex. These are Ward identities relating the tree level vertices and propagators of 



[IV 



the classical Lagrangian. They hold true before including ghost and gauge fixing terms. They are 
still valid at one loop for the pinch technique Green's functions. At tree level, setting D x (k 2 ) = k 2 
in Eq. ( |2.100| ), since \ is massless before quantization, one recovers Eq. Q2.27| ). 



However these Ward identities are in fact irrelevant for the unitarity cancellations that we 
have considered, i.e. those between the W— fusion graphs and the graphs T W f or Tff. They will 
become essential when one attempts to separate out the Higgs graph alone from the rest of the 
W— fusion graphs, and to require for it good high energy behaviour. In this case Dh in Eq. ( |2.100| ) 
will cancel the Higgs propagator and thus extract a piece from the Higgs graph that combines 
with the gauge boson graphs, while the rest of the terms in Eq. (|2.100|) , namely T Hxx and D x , will 
remain in the Higgs graph. This is a cancellation 'internal' to the W— fusion graphs. The extra 
vertex terms in the Higgs graph will however modify its contribution by 0(a). It is conceivable 
that these terms, if not correctly included, may lead again to violation of unitarity. However since 
they are 0(a) this would come about only at extremely high energies. 

In our approach we have neglected all such terms of 0(a) in the numerator of the Higgs graph. 
However, in order to obtain an estimate of the difference between our naive treatment and the 
correct treatment on the resonance, we have calculated the differential cross section using both 
a constant width and the pinch technique running width for the Higgs. We have only found 
negligible numerical differences of O(10~ 3 — 10~ 4 ). These do not affect any of the plots we present 
below. In general we believe that the numerical significance of these terms is very small. 



The squared amplitude of Eq. ( [2.92|) is calculated with FeynCalc and a Fortran output in 



terms of scalar products between momenta and fermion currents is obtained. The phase space 
integration is done by Monte Carlo methods using VEGAS. We decompose the phase space according 
to the structure of the Higgs graph, i.e. to the product of a three body phase space times a two 
body decay of the Higgs: 

G?$ 4 (gi, <?2; <?3, <?4, fa, k 2 ) = (27r) 3 d$ 3 (gi, q 2 ; q 3 , q A , q) dq 2 d$ 2 (q; fa, k 2 ) . (2.101) 



3 The process pp —> ZZ + 2 jets + X 

In this section we investigate the quantitative impact of our results by studying the realistic case 
of heavy Higgs production at the LHC. In particular, we focus on the ZZ final state (the 'gold- 
plated' decay channel for a heavy Standard Model Higgs, see for example Ref. 0) We require, 



in addition, two forward 'tag' jets [11], so that the leading-order subprocess is qq — > qqZZ. We 



will compare the cross sections obtained using the full scattering electroweak amplitude for this 
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process with those obtained using the pinch-approximated W— fusion amplitude, as defined 
and calculated in the previous section. 

The full cross section is obtained by summing over all possible parton subprocesses folded with 
the appropriate parton distributions: 

a(pp ^ZZ + 2j + X)=J2 J dx 1 dx 2 f qi / p (x 1 , Q 2 )f Q2 / P (x 2 , Q 2 )a(qiq 2 -> q[q 2 ZZ ) ■ (3- 1 ) 

For the parton distributions we use the default MRST set | JL2|| , with scale choice Q = Mh- As 
our primary aim is to compare cross sections calculated using different scattering amplitudes, we 
fix the parton and parameter choices throughout the study, and impose the same cuts on the 
final-state particles. Because it is physically indistinguishable, we must include also the resonant 
Z— fusion contribution. The Z— fusion amplitude can be obtained straightforwardly from the 
relevant 3 Feynman graphs, since they form a gauge invariant subset. In some of the subprocesses 
one of the two resonant mechanisms appears as Higgs-strahlung S graphs, which are suppressed 
relative to the fusion graphs by an additional factor of s. In such cases we retain only the 
energetically dominant, fusion one. For example, in uu — > ddZZ the Higgs-strahlung Szz graphs 
are suppressed. When they are both relevant, as in du — > udZZ for example, we will neglect 
the interference between them. This is because the momenta of the final state quarks are crossed 
in the Z— fusion graphs relative to the W— fusion quark momenta. Since each amplitude peaks 
in the forward region there is a negligibly small region of overlap in phase space. For the same 
reason, when calculating amplitudes with identical fermions in the final state we will not consider 
the crossed diagrams, including only the direct one without the symmetry factor 1/2. 

We begin by showing, in Fig. |, the total subprocess cross section as a function of the subprocess 
centre-of-mass energy y/s, calculated three different ways. For purposes of illustration we take 
Mh = 500 GeV. The solid line (F) correponds to the six W— fusion graphs calculated in the 
Feynman gauge. As discussed in the Introduction, this exhibits the unitarity-violating behaviour 
a ~ s 2 at high energy. Isolating the Higgs resonance graph alone (dotted line H) again leads to 
unitarity-violating behaviour, but now o ~ s.[] Finally, applying the pinch technique as described 
in the previous section yields acceptable high-energy behaviour of the form o ~ £n(s) (dash-dotted 
line PT). 

The 'bad' high-energy behaviour of the W— fusion graphs completely swamps the Higgs res- 
onance behaviour. This is illustrated in Fig. which shows the ZZ invariant mass distribution 
at two values of the subprocess energy y/s, again for Mh = 500 GeV. For the W— fusion graphs 
contribution, we see the resonance at 500 GeV disappear under the background as y/s increases. 
At high y/s the cross section is spread approximately uniformly over the range (0, y/s). In con- 
trast, applying the pinch technique gives invariant mass distributions in which the resonance is 
clearly visible at both values of y/s considered. Of course the extent to which the Higgs reso- 
nance is visible over the underlying non- resonant diagram contributions depends on Th- Figure |] 
shows the same ZZ (pinch-technique) mass distribution as in Fig. f|, for Mh = 500 GeV and 
y/s = 4 TeV, together with the corresponding distributions for Mh = 750 GeV and 1000 GeV. 
Notice that, as expected, the resonant peak broadens as Mh increases, becoming barely visible 
over the non-resonant contributions at Mh = 1000 GeV. 

4 Notice also the threshold behaviour of this contribution at y/s ~ M#. 
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Our final two figures show distributions for the full (a/s = 14 TeV) proton-proton cross section, 
i.e. with parton distributions folded in as in Eq. (|3.1| ). In these figures we compare the pinch- 
technique (approximate) result with the full all-diagrams calculation of Ref . || . Figure |] shows the 
ZZ invariant mass distributions for Mh = 500 GeV and Mh = 740 GeV. We see that the resonance 
region is indeed very well approximated by the pinch-technique result. The high- mass tail is also 
in good agreement with the full result. Only the low Mzz region below (9(400 GeV) shows any 
significant difference. Here the large number of non-W— fusion graphs in the full calculation leads 
to an excess over the pinch-technique result. We note, however, that this low- mass region would 
in all likelihood be removed by experimental cuts in a realistic analysis. 

It is relevant to ask whether other features of the final state are well approximated by the 
pinch-technique method. As an illustration, we consider the transverse momentum distribution of 
the forward jets accompanying the Z boson pair. Since it may be possible to detect these jets in 
the LHC experiments, it is important that the pinch-technique result gives an accurate description 
of their production properties. Figure f| shows the predictions for the inclusive jet pr distribution 
in the full calculation and the pinch-technique approximation, for Mh = 500 and 740 GeV. A lower 
Mzz cut is imposed (see above). The distributions are indeed very similar, particularly in shape. 
(The small pinch-technique excess can be traced back to the slightly higher Mzz distribution in 
Fig. |3j.) The rapidity distributions of the jets (not shown) are also very similar. This gives us 
confidence that all important kinematic features of the ZZ + 2 jets production process are well 
reproduced by the pinch-technique method. 

4 Conclusions 

In this paper we have used the pinch-technique to define and calculate a VT-fusion subamplitude 
of the full amplitude for the electroweak process qq — > qqZZ. The pinch-technique amplitude 
is based on the M^-fusion subset of diagrams, and includes the Higgs resonance contribution, i.e. 
WW — > H — > ZZ. More importantly, the pinch-technique amplitude squared is gauge invariant 
and the corresponding cross section exhibits good high energy behaviour. In addition, our numerical 
studies show that the pinch-technique gives an excellent approximation to the full calculation for 
such quantities as the ZZ invariant mass distribution, particularly in the region of the resonance, 
and the rapidity and transverse momentum distributions of the jets accompanying the ZZ pair. 
Not surprisingly, the expression for the pinch-technique is much more compact than that for the 
full amplitude, with a corresponding computation time per event which is between one and two 
orders of magnitude faster. 

In summary, we have demonstrated, in an non-trivial example, that a meaningful and well 
defined separation between signal and background can be achieved even when they contribute 
coherently to the same final state. In the particular case of interest the signal turns out to be 
dominant while its numerical computation compared to the full cross section has been expedited 
enormously. We believe, therefore, that the pinch-technique amplitude could be a useful simulation 
and analysis tool for Higgs production via ly-fusion at the LHC. 
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Figure 2: The TV/ graphs. The relevant kinematics are also shown. Crossed graphs are omitted. 
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Figure 3: The Tff graphs. The graphs not shown refer to permutations of the relative positions 
of the Z bosons on the fermion line to which they are attached. 
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Figure 4: The S graphs. There are 20 S g , Sj and S z graphs respectively. One of the three 
Higgs-strahlung graphs is also shown. 
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(a) 



(b) 



(c) 








(d) 



Figure 6: Pinch terms induced by the momentum of a longitudinal Z boson. The sum of (a), (b), 
(c) and the relevant term for the Tff graphs is a pictorial representation of the Ward Identity of 
Eq. ( |2.30|) when the trilinear vertices in the graphs are the full ones. In the sum, all pinch terms 
on the r.h.s. cancel and the remainder equals the graphs T(x), with the external Z replaced by 
its would-be Goldstone boson \. Terms represented by the ellipsis on the r.h.s. of (b) contribute 
to the remainder R\ of Eq. ( j2.77| ). In (d) we show how pinching is induced in the Higgs graph 
according to the Ward identity of Eq. ( |2.27|) . The pinch term on the r.h.s. contributes to the g a p 
term of F 1 in Eq. ( ^831) . 
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Figure 7: Total subprocess cross section for W— fusion. The solid line represents the six W— fusion 
graphs in the Feynman gauge. The dotted line is the Higgs graph alone. The dash-dotted line is 
the pinch technique result. 



28 




Figure 8: Invariant mass distribution of the Z boson pair. The solid line is the pinch-technique 
result while the dotted line is the result of the 6 W— fusion graphs. 
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Figure 9: Invariant mass distribution for the Z boson pair of the pinch-technique result, for 
different values of the Higgs mass and for subprocess scattering energy y/s = 4 TeV. 
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Figure 10: Invariant mass distribution of the Z boson pair for pp — > ZZ + 2 jets +X at the LHC. 
The solid line represents the full calculation in which all electroweak graphs are included. The 
dashed-dotted line is the pinch-technique approximation in which only the resonant W— fusion 
and Z— fusion graphs are retained. 
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Figure 11: Distributions of the jet transverse momentum for pp — > ZZ + 2 jets + X at the 
LHC. The solid line represents the full calculation in which all electroweak graphs are included. 
The dashed-dotted line is our approximation in which only the resonant W— fusion and Z— fusion 
graphs are retained. 
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